Random Berry Phase Magnetoresistance as a 
Probe of Interface Roughness in Si MOSFETs 
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The effect of silicon-oxide interface rougliness on tlie wealc- 
localization magnetoconductance of a silicon MOSFET in a 
magnetic field, tilted with respect to the interface, is studied. 
It is shown that an electron picks up a random Berry's phase 
as it traverses a closed orbit. Effectively, due to roughness, 
the electron sees an uniform field parallel to the interface as 
a random perpendicular field. At zero parallel field the de- 
pendence of the conductance on the perpendicular field has 
a well known form, the weak-localization lineshape. Here the 
effect of applying a fixed parallel field on the lineshape is 
analyzed. Many types of behavior are found including ho- 
mogeneous broadening, inhomogeneous broadening and a re- 
markable regime in which the change in lineshape depends 
only on the magnetic field, the two length scales that char- 
acterize the interface roughness and fundamental constants. 
Good agreement is obtained with experiments that are in the 
homogeneous broadening limit. The implications for using 
weak-localization magnetoconductance as a probe of inter- 
face roughness, as proposed by Wheeler and coworkers, are 
discussed. 

PACS: 



I. INTRODUCTION 

Disorder has a profound effect on electron transport 
at low temperature. The scaling theory of localization 
applies over an enormous domain and is a keystone in 
our understanding of disorder effects ||l|,D. For two- 
dimensional samples that are weakly-disordered and at 
low temperature ( "weak-localization regime" ) the theory 
predicts the precise dependence of the conductance on 
an applied magnetic field ( "weak- localization lineshape" ) 
The exquisite agreement of the predicted lineshape 
with experiment constitutes an important confirmation 
of scaling theory Q . 

Silicon MOSFETs are an important experimental re- 
alization of a two-dimensional electronic system. In a 
MOSFET electrons are confined to the interface between 
layers of oxide and semiconductor. In this paper we 
analyze the effect of interface roughness on the weak- 
localization lineshape in MOSFETs. Although the effects 
are small they are of interest from various points of view: 

i) The deviations from the known lineshape of an ideal 
interface are small but measurable: our results agree well 
with the experiments that stimulated our work §-§. 



ii) The relevant effect of roughness on the electrons can 
be traced to a subtle quantum interference effect: Berry's 
phase 1^,^. There has been considerable theoretical in- 
terest in designing an experiment sensitive to the influ- 
ence of Berry's phase on quantum electron transport [ p^ . 
Here we have identified a prior experimental detection of 
it. 

iii) A remarkable feature of many quantum transport 
phenomena is their universality, in the sense that the 
observed effects are independent of microscopic sample 
properties such as the mean mobility. The effects of 
roughness have this feature. We find a particularly strik- 
ing regime in which the effect is determined entirely by 
the geometric parameters of the interface and the applied 
magnetic field while being independent of both sample 
mobility and temperature (which controls the dephasing 
length) . 

iv) There has been much experimental and theoretical 
work on the motion of electrons in a random magnetic 
field, particularly in the strong field limit JlT|-p7t (moti- 
vated in part by possible relevance to cuprates and the 
quantum Hall system at filling factor 1/2). We find that 
due to roughness electrons see a uniform in-plane mag- 
netic field as a random perpendicular field; hence it may 
be possible to use this system to study electron motion in 
a random magnetic field. Here our analysis is restricted 
to the weak field limit appropriate for calculating the 
lineshape in MOSFETs. Experiments on other realiza- 
tions of the random field problem are briefly discussed in 
Section VI. 

v) Finally this work may have practical implications. 
MOSFETs are the building blocks of modern electron- 
ics. Roughness of the interface between the oxide and 
semiconductor influences the mobility of the device and 
has been correlated with device failure due to dielec- 
tric breakdown of the oxide ||l^,|l^. It is therefore of 
technological interest to characterize the roughness. R. 
G. Wheeler and co-workers have proposed that magne- 
toresistance measurements on a MOSFET in the weak- 
localization regime can be used as a non-destructive 
probe of its interface roughness [^|J|]. They and others 
have carried out experiments that demonstrate the fea- 
sibility of making the needed measurements The 
present work contributes to this program by providing 
the precise relationship between the interface roughness 
parameters and the measured magnetoresistance. 

Figs. 1 and 2 summarize our findings. Following the 
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FIG. 1. MOSFET with short-ranged correlated roughness 
(1/ -C le)- The conductance (in units of 10~^ mhos) is plotted 
as a function of the applied perpendicular field, Bj_/B^. The 
dotted curve is the classic weak-localization lineshape. It cor- 
responds to zero in-plane field, By = 0. The solid curve shows 
the effect of applying a fixed in-plane magnetic field. It corre- 
sponds to Be/B^ = 400 and B^^ /B^ = {y/8/n^^^)y^Ul/L^. 

experiments of Ref. we imagine the device is placed 
in a magnetic field tilted with respect to the plane of 
the interface. The in-plane component of the field, is 
kept fixed and the conductance is plotted as a function of 
the perpendicular component, B±. To understand these 
curves it is useful to recall that there are two important 
length scales that determine the transport properties of 
the sample: Zg, the elastic mean free path of electrons 
(which determines the mobility), and l^, the distance 
over which electrons maintain phase-coherence and are 
thus able to interfere. In the weak-localization regime, 
lip ^ le- Also, atomic force microscope images reveal that 
statistically the rough interface can be characterized by 
two parameters: A — the root-mean-square height fluc- 
tuations and L = the distance over which the fluctuations 
are correlated. 

For _B|| = we obtain the classic weak-localization line- 
shape appropriate for an ideal interface (dotted curves in 
Figs. 1 and 2). The height|^nd the width of the peak are 
controlled by the dephasing length l,p. The solid curves 
show possible lineshapes when the parallel field is turned 
on. Fig. 1 shows that when the roughness is correlated 
over a very short length scale (i <C /e) the effect of the 
parallel field is to decrease the dephasing length; the line- 
shape is otherwise unaltered. Borrowing the terminology 
of magnetic resonance and atomic physics the effect of 
the in-plane field in this limit may be described as ho- 
mogeneous line-broadening. In the opposite limit when 



^The height is measured relative to the Drude conductance 
given by the large-field asymptotic value of the conductance 
corrected for classical magnetoresistance. 




FIG. 2. Same as Fig. 1 for a MOSFET with long-range cor- 
related roughness (L 3> Ze). The dotted curve corresponds to 
B|| = 0; the solid curve, to By AL = (^sc/10. Also Be = 4005^ 
and L = ^0/^/5. The crossing of the two curves is highlighted 
by a plot of their difference (inset). 

the interface fluctuates slowly {L ^ Zg), the effects are 
more dramatic. The deviation from the ideal lineshape 
changes sign as a function of perpendicular field (Fig. 
2, inset). Again, by analogy to magnetic resonance and 
atomic physics, the extreme limit, L ^ ^ le, can 
be interpreted as inhomogeneous broadening (see Sec- 
tion V) . The sign change is then traceable to a necessary 
inflection point in the ideal weak-localization lineshape. 
The intermediate regime ^ L ^ /g is particularly 
interesting although it cannot be so simply interpreted. 
In this limit the deviation from the ideal lineshape has 
a universal form independent of 1^ and l^. For exam- 
ple, the sign change occurs at a value of the perpendic- 
ular field determined by the purely geometric condition 
B±L^ = (1.79...)/i/e. Note that in all cases the devia- 
tion from the ideal lineshape grows with and A (in 
proportion to i?^ and A^ in the experimentally relevant 
regime) . 

Further discussion and a detailed summary of our re- 
sults is given in Section VII. 



II. BORN-OPPENHEIMER ANALYSIS 

As a model of the silicon-oxide interface, for many pur- 
poses it is sufficient to regard the oxide as an impenetra- 
ble hard wall and to assume that the electrons are bound 
to the interface by a uniform perpendicular electric field 
pof . It will become apparent in the sequel that for the 
present purpose it is only necessary to assume that the 
electrons are firmly bound to the interface: it is not nec- 
essary to commit to a specific form of the confinement po- 
tential C/conf . If the axes are chosen so that the interface 
lies in the plane z — 0, for an ideal interface the confine- 
ment potential would depend only on z; but for a rough 
interface it would vary from point to point and hence 
also depend on x and y. The strength of the confining 
potential controls the extent of the wavefunction in the z 
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direction, denoted i. We also allow for the possibility of 
electron scattering by impurities in the semiconductor by 
introducing a potential t/imp- The magnetic field is as- 
sumed to lie in the x-z plane: thus = B\\,By = and 
Bz = B±. It is convenient to work in a Landau gauge 
and to choose = 0,Ay ~ B±x and Az — B^. The 
Schrodinger equation is 

f d AzY f aY 

-I- C/conf (2, X, y)* + Unnpiz, X, y)* = {x , y, z). (1) 

For later convenience V is used to denote the two- 
dimensional gradient in the x-y plane. Similarly A de- 
notes the x-y component of the vector potential. 

For an ideal interface the electronic motion in the plane 
of the interface and in the perpendicular direction sep- 
arate. The motion in the perpendicular z-direction is 
quantized by the confining potential. Provided the tem- 
perature is sufficiently low and the density of electrons 
not too high, only the lowest subband mode in the z- 
direction is populated and the electronic motion is essen- 
tially two-dimensional. For a rough interface the motion 
no longer separates but we find it is an excellent ap- 
proximation to integrate out the motion perpendicular 
to the interface and to obtain an effective Hamiltonian 
for motion in the plane of the interface using the Born- 
Oppenhcimcr method [pT| . 

For simplicity first suppose there is no magnetic field. 
Assume that the motion in the x-y plane is slow; the 
motion perpendicular to it, fast. The first step in the 
Born-Oppenheimer method is to analyze the motion of 
the fast co-ordinates treating the slow co-ordinates as 
parameters. In this case it is necessary to solve 

- 7;—T^<Pn{z; X, y) + Ucont{z, x, y)(j3n{z] X, y) 

= En{x,y)(l)n{z;x,y) (2) 

Here 4'n{z]x,y) denotes the "local subband" wavefunc- 
tion and En{x, y) is the energy of the n**^ subband wave- 
function. Sometimes it will be convenient to write the 
sub-band states using Dirac notation: 4>n{z;x,y) — > 
\n;x,y). Both the subband wavefunction and the sub- 
band energy vary with location in the x-y plane because 
the confinement potential varies from point to point. It is 
assumed that the subband wavefunction is normalized ev- 
erywhere so that (n; x, y\n; x, y) = dz\(f)n{z] x, y)p = 
1. For a given [x, y) Eq. (2) only fixes the subband wave- 
function up to a phase. It is convenient to choose the 
wavefunctions to be real and to vary smoothly with x 
and y. 

According to the Born-Oppenheimer method an ap- 
proximate solution to Eq. (1) is 



Here 4)g denotes the lowest energy subband wavefunction 
and 4'{x,y) is governed by 

+Eg{x, y)ipix, y) + Wg{x, y)i^{x, y) = EiP{x, y). (4) 

Eq. (4) describes the motion of the electrons in the plane 
of the interface after the transverse motion has been in- 
tegrated out. C/imp, the effect of impurities, now depends 
only on x and y. For impurities that are sufficiently far 
from the interface so that their potential does not vary 
significantly over the confinement scale J/imp = C^imp, 
while for short-range impurities J/imp incorporates their 
effect on electrons in the lowest subband. Note that the 
local subband energy, Eg(x,y) appears as a potential in 
the Born-Oppenheimer effective Hamiltonian. The effec- 
tive potential Wg{x, y) is also determined by the solution 
to the fast problem. It is given by a more complicated 
expression 

Wg{x,y) = ^{g;x,y\{W\n;x,y)) ■ {n;x,y\W\g;x,y). 

(5) 

Next suppose that the magnetic field is turned on. The 
fast co-ordinate is now governed by 

-2^[d-z-'h''^^y) 

-\-Uconiiz,x,y)(l)„{z;x,y) = En{x,y)(j)n{z;x,y). (6) 

Let £,n(z',x,y) denote a normalized solution to Eq. (6) 
when the magnetic field is turned off. This solution is 
chosen to be real and to vary smoothly with x and y. 
It is easy to verify that a solution to Eq. (6) with the 
magnetic field is 



(t>n{z\x,y) 



exp {i^B\\yz^ ^„(z;x,?/). 



(7) 



Turning on the magnetic field leaves the subband en- 
ergy, En{x,y) unchanged; but it introduces a non-trivial 
twist in the subband wavefunctions due to the phase 
factor in Eq. (7). Following Berry this can be 

made explicit by defining the geometric vector poten- 
tial Ag{x,y) = i{Ti/ e){g\x,y\S/\g]x,y) and its curl, the 
geometric magnetic field. Straightforward calculation re- 
veals that the cartesian components of A are 



■^{x,y,z) « ^p{x,y)(j)g{z;x,y). 



(3) 



Agxix^y) = 0; 
^gvix,y) = -B\\Zg{x,y) 

and the geometric magnetic field is 

_ d d 
Bg{x,y) = Q^^gy - q:^^9^ 



(8) 



(9) 
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Here Zg{x,y) is the mean z-coordinate for the lowest 
subband wavefunction, 

/oo 
dzz\^g{z;x,y)f. (10) 
-oo 

Now, according to the Born-Oppenheimer method, an 
approximate solution to Eq. (1) is given by Eq. (3). Here 
(j}g denotes the lowest energy subband wavefunction and 
ip{x, y) is governed by 

%^ / e \2 

V-i-[A(a;,y)+Ag(x,?/)] ipix^y) 

+Eg{x, y)'4){x, y) + Wg{x, y)ip{x, y) 

+t/imp(a;, y)^{x, y) = Eip{x, y). (11) 

Eq. (11) is the central result of this section. It reveals 
that after the fast motion perpendicular to the interface 
is integrated out, the electrons essentially move in two 
dimensions under the influence of an effective potential, 
C^imp + -Eg + Wg and an effective perpendicular magnetic 
field, Bj_ + Bg. 

The effective magnetic field is seen to be the compo- 
nent of the applied field that is perpendicular to the sur- 
face defined by the mean z-coordinate of the subband 
wavefunctions, Zg{x,y) provided that the gradients in 
this surface are small. Note that the geometric magnetic 
field, Bg given by Eq. (9) is proportional to the in-plane 
component of the applied field and to the gradient of 
the surface Zg. It vanishes if = or for an ideal 
interface for which the surface Zg would be fiat. The ef- 
fective potentials Eg and Wg also owe their existence to 
the roughness of the interface but are independent of B\\ . 

Note that Eq. (6) only defines the subband wavefunc- 
tions, up to a phase. A specific choice given by Eq. (7) 
is made in the calculation above. This amounts to choos- 
ing a gauge for the geometric vector potential: a different 
choice of phase would transform the geometric vector po- 
tential, Ag, but would leave the geometric magnetic field, 
Bg unchanged. Thus the calculation above uses a spe- 
cific gauge for both the applied and geometric magnetic 
field — the Landau gauge defined before Eq. (1) and the 
gauge given by Eq. (8), respectively. The gauges are cho- 
sen for their convenience but the results are, of course, 
independent of the choice of gauge. A fuller discussion of 
the gauge invariance of the Born-Oppenheimer method 
is given in chapter 3.7 of Ref. ||]. 

Finally we briefiy discuss the applicability of the Born- 
Oppenheimer approximation. In Appendix A it is shown 
that the approximation should work provided £ ^ Xj, 
XlZg < 1 and £V^Zg < 1. Here i is the typical ex- 
tent of the subband wavefunction in the z-direction and 
A/ is the Fermi wavelength of the electrons. The first 
condition is to ensure that the electron density is not 
so high that more than one subband is occupied. The 
other two inequalities are adiabaticity requirements on 



the surface roughness needed to ensure that the motion 
in the plane of the interface can be approximately decou- 
pled from the transverse motion. They stipulate that the 
roughness must vary slowly on a scale determined by the 
strength of the confining potential. 

To check the validity of the adiabaticity conditions we 
can use the atomic force microscope images of the silicon- 
oxide interface presented in Ref. |^ . These images reveal 
that the height fluctuations of the interface follow a Gaus- 
sian distribution parameterized by Aafm , the root mean 
square height fluctuation, and L, the correlation distance. 
The discussion of Appendix A suggests that the surface 
Zg should be similarly distributed with essentially the 
same parameters. Hence the adiabaticity conditions can 
be expressed in terms of the observable parameters as 
Aafm/-^ ^ 1 and (^Aafm)/^^ ^ 1 and are seen to be 
very well satisfied for both samples of Ref. 1^ . 

III. WEAK-LOCALIZATION LINESHAPE 

The purpose of this section is to review the calculation 
of the weak localization lineshape. The calculation is 
reformulated in a way that is suitable for the eventual 
goal of calculating the tilted field magnetoresistance of a 
MOSFET with a rough interface. 

In the previous section it was shown that effectively the 
electrons move in two dimensions under the influence of 
a random potential, Vix^y)^ and a perpendicular mag- 
netic field. The magnetic field is the sum of the applied 
perpendicular field and the geometric field, which is due 
to interface roughness and the applied parallel field. It is 
assumed that the surface Zg{x, y) is a Gaussian random 
surface with zero mean and variance given by 

(Zg(r)Zg(r'))rough = exp (-^^^) ■ (12) 

As noted above, atomic force microscope images of the 
silicon-oxide interface presented in Ref. Q| reveal a Gaus- 
sian random surface, which makes it extremely plausible 
that the surface Zg(x,y) must also be Gaussian. An ar- 
gument to this effect is given in Appendix A. Eqs. (9) and 
(12) thus determine the statistics of the random magnetic 
field, Bg. 

The random potential V is caused by impurities and by 
the roughness of the interface (in the notation of the pre- 
vious section V — Uimp + Eg + Wg). Here we shall assume 
that the random potential is Gaussian white noise with 
zero mean and variance (V{r)V{r'))i„jp = V^^i^^^^ (r — r'). 
Since we assume that the interface roughness is corre- 
lated over a length scale L, strictly the assumption that 
the random potential is Gaussian white noise is justified 
only under special circumstances (if the impurity poten- 
tial is white noise and it dominates interface roughness 
scattering, or if L is much smaller than all relevant length 
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scales and the impurity scattering is either also white 
noise or is dominated by interface roughness scattering), 
but in practice it is reasonable to believe that our re- 
sults will be more broadly applicable since localization 
effects are believed to be insensitive to microscopic de- 
tails of the random potential. We estimate that for both 
samples studied in Ref. |Q impurity scattering dominates 
interface roughness scattering, and for one sample L is 
shorter than all relevant length scales in addition; but 
it is expected that the results should apply under less 
favorable circumstances also. 

It is useful to consider the electron Green function 
which obeys the Schrodinger equation 



= _^(2)(i.„r'). 



Q{v,v'-E) 
(13) 



For the retarded Green function Q^^ the energy E has a 
positive infinitesimal imaginary part; for the advanced, 
C/"^, a negative part. Linear response theory allows us to 
express the conductance in terms of the Green functions 



jdv jdv'Ag[v,v\Ej) ^ A Ag[v\v,E,). 

(14) 



Here /S.Q = — t/ , Ep = Fermi energy and m = 
effective mass of electrons. We have assumed that the 
sample is rectangular (dimensions L x W) and current is 
driven along the length of the sample which is oriented 
along the y-axis. The two-sided derivative in Eq. (14) is 
defined as 



1-^9= 1^9-9^1- 
ay ay ay 



(15) 



We are interested in the disorder averaged conductance. 
To this end it is necessary to average products of Green 
functions over the random potential and the interface 
roughness under the statistical assumptions made above. 

In the two subsections below we review the calcula- 
tion of the conductance at zero magnetic field and in a 
uniform perpendicular field, circumstances under which 
there is no random magnetic field and it is only neces- 
sary to perform an average over the random potential. 
The effect of the random magnetic field is analyzed in 
the following sections. 



A. Zero Field 

The average over the random potential can be calcu- 
lated perturbatively P,|23[| in an expansion in the small 
parameter (EfTe)^^. Here Te, the elastic scattering time 



FIG. 3. Cooperon diagrams. Solid lines denote impurity 
averaged Green functions. The upper lines are retarded Green 
functions; the lower lines, advanced. Dotted lines represent 
scattering from impurities (shown as crosses). 

for electrons, is given by h / 2tt p{E f)VQ . where p{Ef) = 
density of states for spinless electrons. In this approxi- 
mation the average Green function is given by p3| 

^G«(k,i?/)expzk.(r-r'); 



G^{k,Ef)=[Ef 



h i 

2™ 2t„ 



(16) 



Here G{r,r';Ef) — (C/(r, r'; iJ/))inip is the Green func- 
tion averaged over the random potential. 

The average of a product of Green functions does not 
factorize into a product of the averages; it is correlated 
by the underlying random potential. The correlation be- 
tween and responsible for weak localization is 
called the Cooperon, C(r,r'). Semiclassically this cor- 
relation can be understood to arise from the constructive 
interference of classical paths related by time reversal 
symmetry Fig. 3 shows the Feynman diagrams 

that contribute to the Cooperon 0-0]. From these dia- 
grams we see that it obeys the integral equation 

C(r,r') = C(°)(r,r')+ f dr"C^"\r,r")C{r" ,r'). (17) 



Here C<-°\r,r') = {n/27rp{Ef)Te)\{G^{r,r' , Ef))i^p\\ 
Using Eqs. (14) and (16), the weak-localization contribu- 
tion to the conductance, gwL, can be expressed in terms 
of the full Cooperon: 



2 f 
5WL = ~:^j2^'^e / drC(r, r). 



(18) 



Here D = v'jTe/2 is the electron diffusion constant. Cal- 
culation of 5wL therefore reduces to solution of Eq. (17). 

The customary procedure is to obtain eigenfunctions 
of C'"-* which obey 

J dr'C(0)(r, r')gA(r') = \Qx{r). (19) 

The solution to Eq. (17) is then 



C(r,r') = E 



1- A 



Qx{r)Qxir'). 



(20) 
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It turns out that eigenfunctions that vary slowly on the 
scale of le = VfTe dominate the conductance; hence, it 
is only necessary to accurately calculate the slowly vary- 
ing solutions of Eq. (19). Asymptotic evaluation of the 
Green function, Eq. (16), shows that 



1 



2ttI': 



■ exp 



|r — r 



(21) 



provided |r — r'| ^ A/ and kjle ^ 1- C^"^ is therefore 
sharply peaked about r « r' and for the slowly varying 
eigenfunctions the integral in Eq. (19) can be performed 
approximately by expanding Q\(r') in a Taylor's series 
about r' = r. Keeping terms to second order transforms 
the integral equation (19) into a diffusion equation 



DTeV^Qxir) + Qx{r)^XQxir). 



(22) 



It is sufficient to use the solutions of the diffusion Eq. (22) 
to construct the Cooperon, Eq. (20), and to calculate the 
conductance using Eq. (18). This completes the usual 
calculation at zero field. 

We now introduce an alternative formulation more 
suitable for analyzing the effects of roughness. Consider 
the kernel 



i^(o)(r, r') 



1 



27r/? 



exp 



|r — r 



(23) 



By analogy with C*^"^ and C introduce K given by 
K{r,r') = K'^°\r,r') + j dr" K{y,y")K^^\y" ,y'). (24) 

Note that if^") is not equal to C^") [compare Eqs. (21) 
and (23)]; but the slowly varying eigenfunctions of K'^^^ 
also obey the diffusion Eq. (22). For this reason K and 
C have the same long distance behavior. Thus, we may 
use K instead of C in Eq. (18) to calculate the weak- 
localization conductance since it is dominated by the long- 
distance behavior. 

The chief virtue of K compared to C is that it can 
be expressed in terms of tractable Gaussian integrals. A 
useful expression for K is obtained from Eq. (24) by 
iteration 

oo 

K(r,r') = ^if("Hr,r'); 



K^"\y,y')= j dYi...dYr,K^"\Y,Yi)K'-"'>{Yi,Y2)... 

...K<^'\Yr,,Y'). (25) 

Evidently K^"^") (r, r') may be interpreted as the probabil- 
ity for a Gaussian random walker to go from r' to r in 
n-\-l steps. 

An important physical ingredient must now be added. 
As noted above, weak-localization is a quantum interfer- 
ence effect. It is therefore damped by electron-electron 



and electron- phonon scattering |2j| . In the conventional 
analysis the damping effect of these processes is included 
phenomenologically by adding t^/t^ to the eigenvalue. A, 
in Eqs. (19), (20) and (22). Here is the dephasing time 
and it is assumed that ^ Te in the weak-localization 
regime. The equivalent procedure in our formulation is 
to replace 



i^(")^x(")exp (-^^^^) 



Eqs. (18), (25) and (26) then yield 



9WL 



h L2 



J dr^X(")(r, r) exp 



T0 



(26) 



(27) 



In summary, Eqs. (23), (25) and (27) provide the tools 
to calculate the weak-localization conductance. 

Eq. (27) has an appealing form. gwL is expressed as a 
sum over closed random walks, consistent with the phys- 
ical idea that it is due to the interference of closed paths 
and their time-reversed counterparts. The integral over 
A'(") gives the weight of n-bounce paths; the exponential 
factor shows that the contribution of very long paths is 
cut off by dephasing. For a different formulation of weak 
localization in terms of random walks see Ref. ||2^ . 

Let us explicitly calculate gwL- Successive integration 
over intermediate points shows 



X(")(r 



1 



in +1)1; 



■ exp 



{n + l)ll 



(28) 



a well-known result 
but with L 



/s:("' has the same form as K'-^^ 
in agreement with the general 



\/n + Tle 

proposition that the displacement of a random walker 
grows as the square root of the number of steps. Now 



ffWL 



e2 ^ ^ 1 / nre 
— — > - exp 

n. I , ^ — ^ 11. \ 



=2 

h L 
h L 



=2 W 



Te 

Irh 



n L L 



(29) 



a celebrated result [2 |3|. Here we have introduced l^ — 
^JDt^, the dephasing length, and made use of the in- 
equality T0 ^ Tg. The logarithmic divergence of ijwl 
as Irj, ^ oo can be traced to the slow 1/n decay of the 
weight of long paths. In a field theory formulation it in- 
dicates that the theory flows to strong coupling on long 
length scales but on shorter scales it is "asymptotically 
free" . Provided that the cutoff is shorter than the 
localization length, perturbation theory is accurate. 
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B. Uniform Perpendicular Field 

We consider magnetic fields that are weak compared 
to impurity scattering. The precise condition we assume 
is i? <C Be, where = 0sc/27rZg and 0sc = h/2e is the 
superconducting flux quantum. In this hmit the classical 
paths are not modified by the magnetic field. Its main 
effect is to change the phase of the paths (the Aharonov- 
Bohm effect). This has no effect on the classical Drude 
conductance, but the weak-localization correction, which 
is an interference effect, is suppressed in a manner ana- 
lyzed in Ref. ^. 

The starting point for this analysis is to assume that 
the effect of turning on a perpendicular field on the Green 
function is to multiply it by a phase factor: 

Gg^ (r, r'; Ej) « Gg=o(r, r'; Ej) exp (^z^£ dl ■ A 

(30) 

Here Gb=q is the disorder averaged Green function at 
zero field, Eq. (16); Gb±, the Green function with the 
perpendicular magnetic field turned on; and the integral 
in the phase factor is evaluated along a straight line from 
r' to r. The zeroth order Cooperon becomes 



(0) 



2e 



dl- A 



(31) 



Once again we are interested in the eigenfunctions of 
C^l which obey 

J dr'd°lo{r,r')exp(^t^ J^^ dl ■ A^Q)^{r') = XQx{r). 

(32) 

For the slowly varying eigenfunctions Eq. (32) may be 
transformed into a differential equation by expanding 
both the phase factor and Qx{r') about r' = r to sec- 
ond order. The result is 



DTe V 



i— A(r, 



1 Q,(r)= (A + ^jQA(r). 

(33) 



In Eq. (33) we have explicitly inserted a damping term 
Te/Tcj,. The customary solution is to construct the 
Cooperon from the Landau- level solutions to Eq. (33). 

As in the previous subsection we now describe an al- 
ternative formulation. Consider the kernel 



Kb^ (r, r') = K<^^1 (r, r') + J dr" K^^ (r, r")/^B, (r", r'), 

(35) 

then Kb^^ has the same long-distance behavior as Cs^ 
and may be used in its place to calculate the conductance. 

A useful expression for Kb^^ is obtained from Eq. (35) 
by iteration 

oo 

KBAry)^Y.KB!iry); 

n=0 

K^s'liv, r') ^ j dv,... dr„x|'lo(r, n) . . . 

• ■ ■ ^Eo(r«, r') exp [i^-^ £ dl • a) . (36) 

The integral in the phase factor is evaluated along 
the path from r' to r obtained by joining the points 
r'. r„, . . . , ri, r with line segments in the given order. For 
r = r' 



KBloi^n,r) exp ( «27r-p- 



(37) 



Here = flux through the polygon defined by the closed 
path from r to r and (j)sc — h/2e is the superconducting 
flux quantum. 

Using Eq. (18) and (36) 



5WL — 



h L2 



/OO 
^b!_ (i"' I") exp 



T4, 



(38) 



K'^liry) = 41o(r,r')exp (^^ £ dlA^ (34) 

with iir^^Q given by Eq. (23). The slowly varying eigen- (r,r) 
functions of Kg"^ also obey Eq. (33); hence if we deflne 
Kbj^ via 



In summary Eqs. (23), (37) and (38) provide the tools 
needed to calculate the perpendicular field magnetocon- 
ductance. 

Let us now explicitly calculate the magnetoconduc- 
tance. It is necessary to first evaluate K^^{r,r). The 
integrals are performed in Appendix B and will be used 
in later sections. For the present, more insight is gained 
by noticing that K^j^^ (r, r) has a natural interpretation 
in terms of random walks. 

Imagine that n-sided polygons are drawn on a plane 
starting from a fixed point r. Let the probability density 
of drawing a polygon with vertices at r, ri, r2, . . . , r„, r 
be given by the Gaussian expression [27rZg(n -I- 
l)]x|j°lo(r,ri)i^|3°lo(ri,r2)...i^lj°lo(r„,r). The factor 
2TTll{n+ 1) is included to normalize the distribution. Let 
a denote the directed area of the polygon (defined in Ap- 
pendix C). Then Eq. (37) may be rewritten as 



^-^^^-^yjaF„^.,(a)exp [2m-^ ) . 

(39) 
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Here P„ (a) is the probability distribution of the directed 
area of an n-sided polygon drawn from the Gaussian 
polygon ensemble defined above. Eq. (39) shows that 
i^rj^^ (r, r) is the Fourier transform of the distribution 
of directed areas. The directed area distribution is dis- 
cussed in Appendix C. For n ^ 1 



Pnia) 



and its Fourier transform is 



<^(r,r) 



sinh 



(40) 



(41) 



Eq. (41) is derived in Appendix B and is shown to be 
sufficiently accurate for our purposes for all n provided 
B±^ <C Be- For further discussion of the directed area 
distribution see Appendix C. 

Substitution of Eq. (41) in Eq. (38) leads to 



ffWL 



' h L ^ 2Be 

n=l 



sinh 



nB±_ 
2S7 



exp 



T4, 



(42) 



Since B±^ <^ B^, the sum may be evaluated asymptoti- 
cally.0 Using an integral representation of the digamma 
function [E7| 



we obtain 
ffWL - 



-1 p-xq 

dq[- ^ 



1 W 
2^¥T 



In 



Be_ 

Bi_ 



1 + :^ 

2 Bj_ 



(43) 



(44) 



Here B^ is the field corresponding to a flux quantum 
through the phase-breaking length, B^ = (psc/^T^l^- 
This is the ideal weak-localization lineshape derived by 
Hikami, Larkin and Nagaoka |3| and shown as the dotted 
curves in Figs. 1 and 2. 

It is instructive to summarize the main points of the 
derivation. First the conductance is expressed as a sum of 



The lower limit of the sum requires delicate handling. Intro- 
duce a cutoff > 1 but with NB±/Be < 1 and Nre/r^ < 1. 
The sum over n from 1 to A'' is approximately In A'^ since the 
summand is approximately 1/n over this range. The sum 
from A'^ to infinity may be turned into an integral in the limit 
B±/2Be — > 0. The integral diverges logarithmically at its 
lower limit. The divergence may be extracted by integration 
by parts. This cancels the log dependence of the first sum 
on the arbitrary cutoff A'^ and leaves a finite expression that 
can be expressed in terms of the digamma function using Eq. 
(43). 



closed random walks, Eq. (38). Eq. (37) shows that the 
contribution of a fixed random walk oscillates with the 
enclosed flux with a period (j^scj revealing the Aharonov- 
Bohm origin of the weak-localization magnetoresistance. 
The contribution of all n-bounce paths is the Fourier 
transform of their directed area distribution; it decays 
exponentially for B± ^ (j)sc/nll [see Eq. (42)]. Sum- 
ming over n leads to the gentler decay described by Eq. 
(44). 

IV. SHORT RANGE ROUGHNESS 

In this section we assume that the distance over which 
the roughness is correlated is short compared to the elas- 
tic mean free path (L <C le)- Under this circumstance the 
effect of turning on an in-plane magnetic field is compar- 
atively straightforward: it enhances the dephasing rate 
and hence produces "homogeneous broadening" of the 
weak-localization lineshape. 

For simplicity we first assume that the applied field lies 
in the plane of the interface. According to the analysis 
of Section II, the interface electrons see a uniform in- 
plane field as a random perpendicular field. The typical 
strength of the random field is Bg ~ B^^A/L [see Eqs. 
(5) and (12)]. We assume that the random magnetic field 
is weak compared to impurity scattering (Bg <C Be), a 
condition that is well satisfied in practice. As in the 
previous section, the only effect of the magnetic field in 
this limit is to multiply the Green function by a phase 
factor, so the zeroth-order Cooperon becomes 

4°|(r,r') = 4°lo(r,r')exp (^^^ £ dl.Ag^ . (45) 

The Cooperon is determined in principle by Eqs. (17) 
and (45); our objective is to average it over the interface 
roughness. 

To this end it is convenient to formally solve Eq. (17) 
by iteration: 



Cs„(r,r) = ^Ci7(r,r); 



4°lo(r„,r)exp(z2^'^^ 
dl-A„. 



(46) 



The integral in Berry's phase factor, (f>g, is to be evaluated 
around the perimeter of the polygon obtained by joining 
r, ri, . . . , r„, r in that order. We average the Berry phase 
factor using the gauge given by Eq. (8) and assuming 
that 2g{x, y) is a Gaussian random surface with statistics 
given by Eq. (12). To avoid digression we shall return 
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to the evaluation of this average below. For the moment 

we give the result for L^l^: 



exp j27r— =^ 



rough 



exp I -2^j^BlA^L 



{y-yi? 

|r — ri| 



+ ...+ 



{Vn - yf 



(47) 



From Eqs. (46) and (47) we obtain the average 
Cooperon, 

oo 

CB„(r,r) = 5]4';)(r,r); (48) 

n = 

C^Sl{^,v)= j dri...rfr„4°;(r,ri)...4°)(r„,r); 

4°?(r,r') = cE„(r,r')exp (-2V^^S^A^Lfc^) . 

For brevity {Cb^^ (r, r))rough is written as Cs,, (r, r). 

Evidently Eq. (48) is the formal iterative solution to 
the integral equation 

CB„(r,r')=4°J(r,r')+ j dr"4°)(r, r")CB„ (r", r'). 

(49) 

In principle the roughness averaged Cooperon is deter- 
mined by Eq. (49) with C^°^(r,r') given by Eqs. (21) 

and (48). The expression for C^^(r, r') shows that after 

averaging over roughness, the effect of the in-plane field is 
to multiply the zeroth-order Cooperon by an exponential 
factor. 

To solve Eq. (49) it is most convenient to look for 
eigenfunctions of ■* which obey 



dr'4°lo(r,r')exp 



(y - y' f 



>A(r') 
(50) 



Cb|| may be expanded in terms of these eigenfunctions 
as in Eq. (20). We are interested in eigenfunctions that 
vary slowly compared to le- For these eigenfunctions Eq. 
(50) may be transformed into a diffusion equation by 
expanding Q\ and the exponential factor in powers of 
r — r'. Keeping the leading terms (second order in Q;^ 
and first-order in the exponential factor) we obtain 

(1 + £>reV2) Qxiv) = f A + - + ^"j Qx(v). (51) 

In Eq. (51) wc have put a dcphasing term Tc/t^, by hand. 
The in-plane field is responsible for the Te/ry term. Here 
T]| is defined by 



Til 



.,dy~y') 



l\2 



r — r 



(52) 



Eq. (51) shows that, upon averaging, the only effect of 

an in-plane magnetic field on the long-distance behavior 
of the Cooperon, and hence on the conductance, is to 
enhance the dephasing rate by 



^ -I- 

T4, rii 



(53) 



where ry is given by Eq. (52). 

The form of T|j"^ in Eq. (52) can be obtained from the 
following simple argument. We wish to find the typical 
phase associated with an n-sided random polygon due to 
the random vector potential. Consider a single side of the 
polygon: it has typical length Ig, and the vector poten- 
tial is approximately constant over segments of length 
L but of either sign — A ~ ±i3||A. Then the typical 
phase accumulated on this side is rms[(2e/??,) J A - dV\ ^ 
{2e/h)\Jle/ L)B\\IS.L. Since the phase along each side is 
uncorrelated, the typical total phase is rms[(;/)n-gon] ^ 
^/n{2e/h)^/QL)B\\AL. The contribution of n-gons will 
be cutoff when this phase is of order 1, giving the condi- 
tion for T|j~^. This gives the condition for t^^ if we set 

n T\\/Te. We find Te/T|| ~ 4(e/?i)^B|A^i/e, the same 
form as Eq. (52). 

Now we indicate how to generalize this analysis when 
the applied magnetic field has a non-zero perpendicular 
component. Assuming this component is weak compared 
to impurity scattering {B±^ <C B^), an Aharonov-Bohm 
phase factor must be included in Eqs. (45) and (46). 
After averaging over interface roughness Cb\\,Bj_ obeys 
Eq. (49) with 

= cEo(r, r') exp (i| £ ^I.a) (54) 



exp -2^^BlA?L 



" |r-r'! 



Keeping track of the phase factor, we find that the slowly 
varying eigenfunctions of C^'' obey Eq. (33) with the 
substitution given in Eq. (53). 

Hence we arrive at the central result of this sec- 
tion. For short-range correlated roughness, application 
of an in-plane field increases the dephasing rate accord- 
ing to Eq. (53) hut does not otherwise alter the weak- 
localization lineshape. The dependence of the conduc- 
tance on the perpendicular field is still given by Eq. (44) 
but with 



Ba, ^ Bfh -\- TT^/^ 



h L. 



(55) 
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Since the form of the hneshape is not changed, this ef- 
fect is analogous to homogeneous broadening in magnetic 
resonance and atomic physics. 

The idea that an in-plane field produces homogeneous 
broadening of the line was anticipated by the authors of 
Refs. Our result differs in two important respects: 

First, we find that homogeneous broadening occurs only 
if the roughness is correlated over short distances (L ^ 
le)', and second, the dependence of T|| on the interface 
roughness parameters that we derive [Eq. (52)] is quite 
different from that conjectured in Ref. 

Fig. 6 shows that Eqs. (44) and (55) provide a good fit 
to data with essentially no adjustable parameters. More 
details are given in Section VI. 

To complete the analysis we must now average Berry's 
phase factor exp(i27r(/)g/(/)sc) over interface roughness to 
obtain Eq. (47) . Recall that the distribution of a linear 
combination of correlated Gaussian random variables is 
also Gaussian. The distribution of 4>g is therefore Gaus- 
sian. Evidently it has zero mean since we have assumed 
{Zg{x,y)) rough — 0. For a Gaussian random variable, s, 
with zero mean it is easy to show 



(exp(zes)) =exp(-ie'(s'))- 



(56) 



Hence averaging Berry's phase factor reduces to the cal- 
culation of the variance of the phase. Explicitly we must 
calculate 



^g/roug 



n+1 71+1 „ „ 

^ = T.T.{jdh-A,Jdl,- A,)^^^^^. (57) 

2—1 j — 1 



Here J dli ■ Ag denotes the line integral along the edge 
joining r^-i to (with the understanding that ro = 
r„+i = r). 

Let us first evaluate the diagonal terms in Eq. (57). 
Let the i^^ edge make an angle 9 with the x-axis and let 
s be the length parameter along that edge. Using Eq. 
(8) for Ag and Eq. (12) for the correlation of Zg{x,y) 
we obtain 



( j dh-Ag)'' 



rough 

B?A2sin2 6' / ds j ds'exp 

Jo Jo 

2A2r - y^-lf 



\n - ri_i| 



(58) 



We have assumed that the edge is much longer than the 
correlation length L. Since the length of a typical edge 
is le, this is justified in the short range limit L <^ le- We 
have also used jr^ — ri_i| sin 6 = yi — Vi-i- 

Evidently the only important off diagonal terms in Eq. 
(57) are those for which the two edges pass within a cor- 
relation length of each other. The contribution of such 
edge pairs to the sum will be smaller than the diagonal 



terms, typically by a factor of L/le (except in the rare 
circumstance that the two edges meet at a glancing angle 
and overlap for a considerable part of their length). Fur- 
thermore, the number of contributing off-diagonal terms 
is comparable to n, the number of diagonal terms; hence 
the contribution of the off-diagonal terms may be safely 
neglected in the short range limit L ^ l^- 

That the number of important off-diagonal terms is 
comparable to n can be established by an argument com- 
monly used in polymer physics [ p8[ . Suppose that line 
segments of typical length Ig are drawn at random on a 
plane with a density p segments per unit area. If a fresh 
segment of length 1^ is now drawn on this plane, typi- 
cally it will suffer pl^ intersections. Working in units 
where 1^ — I, the area of a typical polygon in Eq. (46) 
is n and hence the density of edges p ~ n/n = 1. Each 
edge therefore typically encounters one intersection and 
the total number of intersections is comparable to n. 

Adding the diagonal contributions [Eq. (58)] and ne- 
glecting the off-diagonal ones we obtain 



{yi - yf 

Iri -r| 



{y~ynf 



(59) 



Use of Eqs. (56) and (59) finally leads to Eq. (47). 



V. LONG RANGE ROUGHNESS 

We now consider the roughness correlation length to 
be long compared to the mean free path [L ^ /g). This 
circumstance is more difficult to analyze; hence, the dis- 
cussion is limited to the circumstance that the in-plane 
magnetic field is weak in a sense made precise below. 
Roughly it is assumed that the random Berry phase is 
small for all contributing paths. ^ The experiments of 
Refs. meet this condition for most part. For simplic- 
ity, we first assume the applied field lies in the plane of the 
interface, B+ ~ (Section Va). In the next subsection 
no restriction is placed on B+ (other than B+ <C Be). 



A. Zero Perpendicular Field Applied 

Since Bg <Si Be is automatically ensured by the weak 
field condition on that we will impose below, the 
Cooperon is given by Eq. (45). Since the Berry phase 
factor varies on the scale of L ^ Zg, we may work with 
K instead of C and write 



■^This is more stringent than the condition Bg ~ A/L ^ 
Be imposed in the previous section [see discussion following 
Eq. (45)]. 
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K^^l (r, r') = Kf^,{v, v') exp (^z| £ d\.A, 
in place of Eq. (45) and 

oo 

^i",^ (r, r) = y rfn . . . rfr„ii'^°lo(r, n) . . . 

• • • i^slolrn, r) exp ( «27r^ 

V ^sc 



(60) 



(61) 



in place of Eq. (46). Note that we could not make this 
replacement in the previous section since the Berry phase 
fluctuated rapidly on the scale of le- Our objective now 
is to average Kb^^ over interface roughness and then use 
Eq. (27) to evaluate the conductance. 

The Berry phase factor is a Gaussian random variable 
so it is easy to perform the average. Using Eq. (8) for 
the vector potential, Eq. (12) for its statistics and Eq. 
(56), we obtain 

/C(;;)(r,r) = j dr^...dr„K^^Ur,r,)... (62) 
• • • K]sLo{j^n, r) exp [-V{r, ri, . . . , r„)] 



with 



y(r,ri,...,r„) = 

■ z^ivj - yj--^)iyk - Vk-i) 

i=l fe=l 



(63) 



X exp 



-^(ri +rj-i -rfc -rfc_i) 



Here /C^") = rough and tq = r„+i = r. /C^") is 

thus an integral over (n+ l)-bounce closed random walks 
just as it is at zero field but the weight of each polygon 
is no longer Gaussian. The links of the polygon now 
"interact", and the interaction, given by Eq. (63), is 
anisotropic and long-ranged. 

/C^"^ is difficult to calculate. For weak in-plane mag- 
netic fields it is sufficient to analyze the interaction per- 
turbatively. To first order 



<54';;(r,r)^4';)(r,r)-<io(r,r) 

= - j dri... drnK^eloi^' ri) ■ • ■ 



(64) 



. . . K^eloi'^n, r)V{r, ri, . . . , r„). 

V{r, ri, . . . , r„) is a sum ovcir pairs of links according to 
Eq. (63). To illustrate the evaluation of Eq. (64) focus on 
the particular term for which j = 1 and A; = 3, 4, . . . , n. 



It is convenient to perform the integrals over the end 
points of the j"^ and fc"^ links last. Integrating over the 
remaining intermediate points yields 

(r, r) = j drdridrk-idrk (65) 

X ^ z — \-rT-. — —{yi - y){yk - Vk-i) 

(rfe_i - ri)2 



(n + 1- A:)(A:-2) 
(r-rfe)2 



X exp 
X exp 
X exp 



2ll{n + l-k) 2ll{k - 2) 



(rfe - rfe_i) 



(ri 



4L2 



2/2 



(r -I- ri - rfe- 



This is a low-dimensional Gaussian integral and can be 

explicitly evaluated. 

Proceeding in this manner we obtain 

<54^,Hr,r) = + /2(n,d) + /3(n,d)]; 

n(l + nd) 



fi{n,d) 



/2(n,d) 
/3(n,rf)«l 



16rf2 l^^[k'^ -nk-{n/Ad\Y' 

\+nd 



(66) 



with d = /g/(2L2). The three contributions to /C^") arise 
from the interaction of disconnected links, adjacent links, 
and from the self- interaction of links, respectively. We 
have used 1 and d <C 1 to simplify Eq. (66). This is 
justified because weak-localization is dominated by long 
paths and we arc concerned with long-range correlated 
roughness in this section. 

Using Eq. (27) for the conductance and Eq. (66) for 
, after considerable simplification we obtain the par- 
allel field magnetoconductance 



^9 = -'^UB,ALfm 



with a = li/L^ and 



/(a) = a / 
Jo 



dxe 
X I 1 + 



(67) 



(68) 



1 ^ ^/T+xa — ^/xa\ 

2[xa{l + xa)]^/'^ ^l + xa + ^/xa J ' 



It is striking that Sg is independent of the mean free path 
le- In appropriate units it is a product of two factors: 
the flux through an area AL, determined entirely by the 
geometry of the rough interface, and f{a). f{a) describes 
the crossover as is varied relative to L. It is plotted in 
Fig. 4 and has the asymptotic behavior 
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f(a) / a 




FIG. 4. Plot of /(a)/a as a function of a = ll/L'^- /(«) de- 
scribes the dependence of the in-plane magnetoresistance, at 
zero perpendicular field, on the roughness correlation length, 
L (see Eq. 67). 

f{a) w a for a 3> 1; 

2 2 , 
~ —a LOT a <C 1. 

(69) 

Substituting Eq. (69) in Eq. (67) gives the asymptotic 
behavior of the parallel field magnetoconductance 



2;2 



55: 



e2 B'A'l^ 



3 n 



for > L; 



for L > /rf 



(70) 



These formulae have a simple interpretation. The ran- 
dom magnetic field has a typical magnitude A/i and 
is correlated over a distance L. For L l^, the typi- 
cal flux through an area is therefore (ByA/L)?^; for 
L <C it is AZ^ (see footnote^). Eq. (70) then shows 
that Sg (in units of e'^/h) is the square of the flux through 
a phase-coherent region in units of the flux quantum. 



B. Perpendicular Field Applied 

The calculation with a perpendicular field applied 
closely parallels the calculation at B± = 0. We write 



K 



(0) 

B I ,B\\ 



.2e 



£dl-{Ag + A) 



(71) 



in place of Eq. (60) and 



*The area can be broken up into l^/L^ correlated 
squares of area L^. The typical flux through each square 
is (_B|| A/L)L^. Since it can be of either sign, the flux through 
the total area grows as the square root of the number of 
squares. 



K 



B±,B, 



(r,r) 



4n) 



00 

n=0 



(r,r); 



K^bLb, (r, r) = 1 . . . dr„ifi°lo(r, n) 



exp ^i27r^-^ ) exp ( j27r^^ 
dl • Ag; 
dl-A 



(72) 



in place of Eq. (61). Here (j)± is the Aharonov-Bohm flux 
through the polygon that goes from r to r via ri . . . r„ 
due to the applied perpendicular field; (pg is the Berry 
phase around the same polygon. 

The Berry phase factor is a Gaussian random vari- 
able. Performing the average as in the previous subsec- 
tion yields 

I^^bIb, (r, r) = 1 dn . . . dr„x(,°lo(r, n) . . . 

■ • • ^sio('^"' i") exp[-y(r, ri, . . . , r„)] 



X exp i27r 



(73) 



Here /C(") = {K'^'''>) rough and V is given by Eq. (63). 
As in Section lllB, /C'") is the Fourier transform of the 
directed area distribution of n-sided polygons; but the 
weight of the polygons is no longer Gaussian. The links 
of the polygons now "interact" and the interaction given 
by Eq. (63) is anisotropic and long-ranged. 

For weak in-plane fields a full analysis of /C*^"-' is not 
needed. It is sufficient to analyze the interaction pertur- 
batively. To first order 



(n) 

B±,B„ 



r,r)^4"^B„(r,r)-n"/(r,r) 



■'B±,B 



j dvi...dVr,K^gl^{v,Vi) 



.(0) 

X exp I 



^B=o(i-n,r)T/(r,...,r„ 



(74) 



V{r, ri, . . . , r„) is a sum over pairs of links according to 
Eq. (63). To evaluate Eq. (74) it is convenient to focus 
on the particular term corresponding to the j*'' and k^^ 
links. As in the corresponding evaluation of Eq. (64), 
it is convenient to perform the integrals over the end 
points of the j^^ and k^^ links last. The integral over 
the other intermediate points is facilitated by Eqs. (BIO, 
B15) of Appendix B. The expression that results is a 
low-dimensional Gaussian integral over the end points of 
links j and k that can be explicitly computed. 
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Proceeding in this manner we obtain 



eBuA 



h 



[/1 + /2-/3]; 



/i(n, 7, « — Xl{cotli[fcV7] coth[(n - k)^ - 1 

fe=3 

+ coth[fc^ + coth[(n - k)^/y\} 

X {sinh[fc^ sinh[(n - k)^/y]}~^ 

X (coth[fcV7] + coth[(n - fc)y^] + 4-^) 



/2(n,7,rf) 



2n V7 
TT sinh[n^/7] 



{^coth[nv^ + rf} 



?T 1 

fain, J, d) « -Vr^-Tr — FT 



(75) 



where 7 = {B±/2Be)'^ is a measure of the dynamical 
strength of the appUed perpendicular field and we have 
tafcen n 1, 7 <C 1 and d <C 1. This is justified because 
weak-localization is dominated by long paths, the applied 
field is dynamically weak, and we are concerned with 
long-ranged roughness in this section. We have also taken 
n7 <C 1, justifiable a posteriori, because we find SK.^"'^ 
decays exponentially for n ^ 1 /y^. 

Using Eq. (27) for the conductance and Eq. (75) for 
^^B_^,S|i' after considerable simplification, we obtain a 
lengthy expression for the parallel field magnetoconduc- 
tance. 



dg = -2 



e2 /BuAL 



Here rj 



/2-kB^L'^ and 



u{t,,B^/B^). 



(76) 



u{r],B^/Bj_) = 



dy- 



Iny 



y 



v{y,r]) 



(77) 



with 

v{y,r]) 



2 [(l + 27? + 2?7^)(y^-l)-H4r?] 



[(1 + 2r?)(y2 _ 1)2 + (8 + 4r]){y^ - 1) + 8] 



p/2 



f [^/7-2r?(y2-l)]2-(^2_l)2 

I[y7+2r?(y2-l)]2_(2;2_i)2 

f{y,v) = iy' - 1)[(1 + 2r?)'(y' - 1) + 87?]. 



(78) 



Eqs. (76)- (78) are the main results of this section. 
They give the shift in conductance from the standard 
weak-localization lineshape due to the applied in-plane 
field Bu. Remarkably Sg is independent of le- 




10 



h (T1) 

FIG. 5. Plot of h{'n)/r] as a function of ?? = <f)sc/{2nL^B±). 
The inset shows the small r) behavior and the zero of h{r]) 
in detail. h{r]) describes the perpendicular field dependence 
of the in-plane magnetoresistance in the geometric limit, 
le«L«l^; see Eq. (79). 

Since the expression for 6g is complicated, it is instruc- 
tive to examine several special cases. First, by taking the 
limit B± ^ it is possible to recover the result of the 
previous subsection, Eq. (67). Next consider the strong 
perpendicular field limit B±E ^ ^sc- In this limit 



e2 /B||AL\ 
^9 = -2- ( hiv) 



n 



Psc J 



with 



Hv)=vJ dylny v{y,r]) - —rj. 



(79) 



(80) 



A plot of ft(ry) is shown in Fig. 5. It has the asymptotic 
behavior 



h{r]) Pi T] for ?7 <C 1 
« — — for ri~:^ 1. 



(81) 



Note that /i(r/) changes sign; it vanishes for « 0.558. 

Eq. (79) reveals that for B±^P^ ^ (f)sc, Sg depends only 
on the magnetic fields and B± and geometric param- 
eters that characterize the rough interface, A and L. It 
is independent not only of l^ but also ofl^. The interme- 
diate regime in which l^, <C L, but L <C l,j,, is particularly 
interesting. In this limit B^ <C (psc/L^ and hence Eq. 
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(79) applies in a range of B± over which t] varies from 
large to small. Thus the perpendicular field magnetocon- 
ductance curves at different values of By cross at a point 
determined by the purely geometric condition 



B±L^ w 1.79-. 



h 

)-. 

e 



(82) 



For this reason we call this intermediate regime, I 
L <^l^, the geometric regime. 

The opposite limit of extremely long-ranged rough- 
ness, L 3> ^0 ^ Ze, we call the inhomogeneous broad- 
ening regime. In this limit, too, perpendicular magneto- 
conductance curves corresponding to different small 
cross, but the crossing point is not determined by the 
purely geometric condition, Eq. (82). In this case the 
range of B±^ over which the simplified formula, Eq. (79), 
applies corresponds to very small rj and does not include 
the geometric crossing point, rj — 0.558 . . . We must 
therefore use the general expression Eq. (76) to deter- 
mine the crossing point. Nonetheless in the inhomoge- 
neous broadening regime the deviation from the conven- 
tional weak-localization lineshape and the existence of a 
crossing point has a simple interpretation. Since the ran- 
dom magnetic field is correlated over a length scale i, the 
sample breaks up into blocks of size each of which sees 
a slightly different (but uniform) perpendicular magnetic 
field, B±+SB. For small SB, the shift in conductance of a 
particular block may be computed from the conventional 
weak-localization lineshape g-wL{B±) by expanding in a 
Taylor series: 



Sg 



d 1 

gwL{B^)SB + - " 



dBjL 



2dBl 



g^i^{B^)5B^ 



(83) 



Bearing in mind that the average value of the random 
field is zero and the typical value, B^^A/L, the average 
deviation is 



Sg^B; 



J^dBl 



gwL{Bi_). 



(84) 



The curvature of gwL{B±) is positive at B± = 0. At large 
B± > (f>sc/ll, 3wl(S_l) ^ {s'^ /h)hi{B±/ B^); hence the 
curvature is negative. Thus there must be an inflection 
point at which the curvature of g^^i^{Bj_) vanishes. Ac- 
cording to Eq. (84) this is the crossing point. Substitu- 
tion of the known asymptotic forms of gwh{B±) in Eq. 
(84) allows a check on the consistency of this interpreta- 
tion. For small B±^ we recover the estimate given in the 
second line of Eq. (70). For B±^ ^ 4'sc/l'i we obtain 



Sg 



h 52 L2 



(85) 




in agreement with the small 77 limit of Eq. (79). 



B 



FIG. 6. Plot of the conductance (in units of 10"'' mhos) 
against perpendicular field, Bx (in Gauss), for the control 
device of Anderson et al. n. Data points for By = are 
grey; for B\\ = 1.5 T, black. The dashed line is a fit to 
the weak-localization lineshape, from which Anderson et al. 
infer — 0.085/im, = 0.76/xm. The solid curve shows 
the fineshape at _B|| — 1.5 T calculated from the expression 
for homogeneous broadening derived in this paper, Eqs. (44) 
and (55). All the parameters for this curve are independently 
determined. The values of A = 1.35A and L = 0.055/im are 
obtained from atomic force microscope images; for 1^ and l,/,, 
we use the values cited above. 

VI. COMPARISON TO EXPERIMENT AND 
RELATED WORK 

In their experiments Wheeler and collaborators apply 
a fixed -By and measure the conductance as a function of 
B±. They have reported such lineshape measurements 
for many devices and have also explored the effect of 
varying electron density on the lineshape. The data of 
Anderson et al. Q are particularly interesting because 
they have independently measured the interface rough- 
ness by etching away the oxide and imaged the exposed 
silicon surface using atomic force microscopy. Hence we 
focus on the measurements reported in that paper. 

Fig. 6 shows the measured magnetoconductance of the 
control device at JSy = (grey) and By = 1.5 T (black). 
By fitting the By = data to a weak-localization line- 
shape (dashed curve in Fig. 6), Eq. (44), Anderson et al. 
infer l^ = 0.085 //m, l^ = 0.76 /im. Atomic force micro- 
scope measurements yield A — 1.35A and L — 0.055 /im. 
Thus the device is in the short-range correlated regime. 
The black curve shows the calculated lineshape using the 
independently measured values of A, L, 1^ and l^. No pa- 
rameters are adjusted except for a constant offset. The 
fit is good, comparable, for example, to the fit to weak- 
localization lineshape at B\\ = 0. 

Anderson et al. also report measurements on a sec- 
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ond device (called the textured device). This device is 
fabricated on a silicon wafer that has been roughened by 
argon sputtering prior to gate oxidation. As a result the 
textured device has atypical values of A and L. L for the 
textured sample is comparable to le- Consistent with this 
we do not obtain a good zero-parameter fit to either the 
short or long-range roughness expressions derived here. 
Following Anderson et al. we can force-fit the data to 
the homogeneously broadened lineshape; but the value 
of A^L we then obtain is inconsistent with the atomic 
force microscope images. 

Finally it is interesting to compare the experiments 
of Wheeler et al. to other work on random or periodic 
magnetic fields. Such fields have been realized by gat- 
ing a GaAs/AlGaAs heterostructure with type II super- 
conducting film [ PIJT2I ; by randomly depositing type I 
superconducting grains or ferromagnetic dysprosium 
dots 1 14 1 on the heterostructure surface; and by cover- 
ing the device with a rough macroscopic magnet 
In work most closely related to Wheeler and co-workers, 
Gusev et al. have prepatterned the substrate by drilling 
a lattice of holes so that the interface at which the two- 
dimensional electron gas forms is periodically modulated 
p6t . An inhomogeneous magnetic field leads to many 
interesting features in the magnetoresistance even in the 
classical and ballistic regimes; many experiments have 
focused on these effects. Except for the work of Rammer 
and Shelankov |jl^ noted below, much of the theoret- 
ical literature too explores these regimes or focuses on 
the fundamental question of whether all eigenstates are 
localized in a random magnetic field. 

Two experiments study quantum corrections in the 
weak-localization limit: In the work of Bending, von Kl- 
itzing and Ploog [ pT| the inhomogeneous magnetic field 
is produced by an Abrikosov lattice of flux lines in the 
superconducting gate; hence it is periodic. However the 
magnetoresistance is not sensitive to the arrangement of 
the vortices since the dephasing length is shorter than the 
distances between the vortices. This circumstance was 
analyzed by Rammer and Shelankov |r^ . In the experi- 
ment of Gusev et al. too the magnetic field is essentially 
periodic rather than random allowing them to analyze it 
via a chessboard model jl^. Thus in both these works 
the inhomogeneous magnetic field is rather different from 
the random Gaussian correlated field studied in the ex- 
periments of Wheeler and co-workers and analyzed here. 



VII. SUMMARY AND CONCLUSION 

In this paper we have been concerned with the in- 
fluence of interface roughness on the magnetoconduc- 
tance of MOSFETs. Using a Born-Oppenheimer ap- 
proximation we have shown that effectively electrons re- 
spond to an applied in-plane magnetic field as though 
it were a random perpendicular field. Technically, the 



random magnetic field appears as a Berry phase term in 
the Born-Oppenheimer effective Hamiltonian. Eq. (9) 
shows that the random magnetic field is the component 
of normal to the rough surface Zg{x,y). The Born- 
Oppenheimer analysis also gives precise meaning to the 
surface Zg{x,y): it is the location of the centre of the 
local subband wavefunction. 

The next step is to analyze the magnetoconductance. 
This is controlled by the familiar physics of weak- 
localization: at zero field there is a contribution to the 
conductance due to constructive interference of closed 
electron paths and their time-reversed counterparts. This 
contribution is suppressed by application of a perpen- 
dicular magnetic field because the electron acquires an 
Aharonov-Bohm phase as it traverses a closed path. The 
analysis of this paper shows that an in-plane field too 
has an effect due to interface roughness. The electron 
acquires a Berry's phase equal to the flux of the random 
field through the closed path. The effect of the random 
Berry's phase on the conductance is the main focus of 
this paper. 

To calculate the effect it is necessary to sum the stan- 
dard divergent series of Cooperon diagrams (Fig. 3). We 
find it useful to reformulate this sum in a manner remi- 
niscent of the central limit theorem. Feynman graphs of 
high order control the sum. In these graphs we are able, 
without significant error, to replace pairs of Green func- 
tion lines with suitable Gaussian factors. At zero field, 
the 71 order diagram is then essentially the weight of a 
closed rt-step Gaussian random walk [Eq. (27)]. For the 
known case of an applied perpendicular field, this for- 
mulation brings out clearly the relationship between the 
classic weak-localization lineshape (Eq. 44) and the dis- 
tribution of directed area for closed random walks (see 
Eqs. 38 and 39 and Appendix C). The method also 
proves useful when an in-plane field is applied for long- 
range correlated roughness (L ^ ?e). In this case the 
magnetoconductance is related to the directed area dis- 
tribution of an interacting random walk (see Eqs. 63 
and 73). To make contact with experiment it is sufh- 
cient to analyze this problem pcrturbatively. We do not 
attempt a complete solution of the interacting random 
walk or calculation of the full magnetoconductance for 
long-range correlated roughness. These problems are left 
open for future work. 

For short-range correlated roughness, however, we are 
able to obtain the full magnetoconductance. The depen- 
dence of the conductance on perpendicular field, B±, still 
has the classic weak-localization lineshape (Eq. 44) when 
an in-plane field is applied, just as it does at By — 0. The 
only effect of the in-plane field is to enhance the dephas- 
ing rate (Eq. 52). In terms of field scales, becomes 
larger; the effective is given by Eq. (55). Hence we say 
that short-range correlated roughness produces homoge- 
neous broadening of the weak-localization lineshape. 

The departure from the conventional lineshape is more 
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striking in case of long-range correlated roughness (al- 
though quantitatively smaller if all parameters except L 
are held fixed). The change in conductance, 8g, when a 
parallel field is applied ("parallel field magnetoconduc- 
tance") is given by Eq. (76). The most remarkable fea- 
ture of this expression is that it is independent of 1^ ■ Sg 
is determined entirely by the magnetic field, geometric 
parameters of the interface, A and L, and the dephas- 
ing length, l^. It may be of either sign (positive near 
B± ~ and negative for large B±). The expression for 
6g is lengthy and cannot be more compactly expressed in 
terms of known special functions. It can be numerically 
computed and plotted with ease and precision; however 
it is worthwhile to examine several special cases. 

In the limit of extremely long-ranged roughness, L ^ 
lip 3> le, the parallel field magnetoconductance has a 
simple interpretation: the sample breaks up into inde- 
pendent blocks of size each of which sees a slightly 
different magnetic field. Hence in this limit we say that 
the effect of an in-plane field and roughness is to inhomo- 
geneously broaden the weak-localization line. Pursuing 
this interpretation we can understand the sign change in 
6g and recover its form in various circumstances [see the 
discussion following Eqs. (70) and (82)]. 

The limit of intermediate-range correlated roughness, 

^ L ^ le, is particularly striking. For sufhciently 
large perpendicular field, B±l'^ 0sc, Sg becomes inde- 
pendent not only of le but also [see Eq. (79)]. Essen- 
tially the parallel magnetoconductance is controlled by 
the function h{rj) where 77 — (l)sc / {2tt L'^ B ±) is a measure 
of the perpendicular flux (see Fig. 5). 

Another special case, useful for making contact with 
experiment Q] is the parallel field magnetoconductance 
for B± = (Eq. 67). In this case, Sg, in units of e^/h, 
is a product of the squared flux (i3||A_L)^, determined 
entirely by the field and the geometry of the interface, 
and the function f{a). f{a) describes the crossover as 
is varied relative to L {a = I'^/L? . See Fig. 4). 

Wheeler and co-workers have shown experimentally in 
a number of cases that the effect of an in-plane field is to 
produce homogeneous broadening. Our analysis of the 
short-range correlated regime provides a theoretical jus- 
tification of such fits to data and allows quantitative in- 
formation about the interface roughness to be extracted. 
It is encouraging that for the control sample studied by 
Anderson et al. [^, for which the interface parameters 
are independently measured via atomic force microscopy, 
our formulae give a satisfactory zero parameter fit (Fig. 
6). Unfortunately, our analysis shows that in the short- 
range regime, magnetoresistance measurements will not 
separately yield A and L: instead they provide the com- 
bination A^L. Thus it will be necessary to combine mag- 
netoconductance measurements with other techniques to 
separately obtain A and L. 

The roughness of the interface is controlled by its pro- 
cessing. For example, Anderson et al. [o] are able to 



increase both A and L by argon sputtering the silicon 
surface prior to gate oxidation. It would be desirable 
to create devices with long-range correlated roughness in 
order to experimentally observe, for example, the geo- 
metric regime described above. At the same time, for 
long-range correlated roughness, the prospects are much 
better for extracting both A and L from magnetoconduc- 
tance measurements alone, provided the measurements 
can be made with sufficient precision. Eq. (76) provides 
the means to study the feasibility of such experiments. 

Another potential application of our analysis is related 
to efforts to engineer the MOSFET interface to have spe- 
cific structure, for example, a structure that is periodic 
in one direction like a corrugated sheet ||2|]. Weak- 
localization magnetoconductance measurements would 
provide a non-destructive way to check that structures 
have been successfully fabricated. It should be a straight- 
forward and interesting extension of our analysis to cal- 
culate the magnetoresistance signature of various simple 
periodic structures. 
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APPENDIX A: ADIABATIC APPROXIMATION 

Here we briefly discuss the conditions for the validity of 
the Born-Oppenheimer approximation used in Section II. 
Roughly, we want only the low-lying states of the lowest 
subband to be occupied and we want the gap between 
subbands to be large. For this purpose, and to ensure the 
two dimensionality of the electron gas, we need the Fermi 
energy to be small compared to the subband spacing. 
The subband spacing is of order 'h'' jml? where t is the 
length scale over which the electron is confined. Hence 
this condition may be phrased as 

A/ » L (Al) 

To obtain additional conditions we return to the 
Schrodinger Eq. (1) and seek a solution of the form 

= ^^/'n(a;,?;)(?;'„(2;;a;,y). (A2) 

n 

For simplicity we limit discussion to the case of zero mag- 
netic field, ipnix^y) then obeys Pj2l|] 
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-—V'^ + Uimp{x,y) + En{x,y) + Wn{x,y) 
2m 



i^n{x,y) 



(A3) 



Here 



n h 
Hnm = — {n;x,y\W\m;x,y).W + — {n; x , yl^"^ \n; x , y) . 
m Zm 



Eq. (A4) constitutes an exact reformulation of the 
Schrodinger Eq. (1). The Born-Oppenhcimer approx- 
imation consists of keeping only one term in the sum 
(A2) and hence omitting the off-diagonal terms in (A3). 

In a commonly used model of the interface |20| it is 
assumed that the confining potential has the form 



Uconiiz; X, y) = U[z - C,{x, y)]. 



(A5) 



Here C,{x^ y) is the elevation of the silicon-oxide bound- 
ary. In this model the mean elevation of the local sub- 
band wavefunction exactly tracks the interface; hence the 
surface Zg{x, y), probed by weak-localization magnetore- 
sistance, has the same statistics as the silicon-oxide in- 
terface ^{x,y), which is presumably imaged by atomic 
force microscopy. For a more general form of the con- 
finement potential this precise correspondence between 
Zg{x,y) and Ci^iU) is lost but it seems reasonable that 
the two surfaces would have similar statistics, compara- 
ble correlation lengths and mean-square fluctuations. 

Eq. (A5) implies that the local subband wave function 
is of the form 



(j)n{z]x,y) = ^n[z - C{x,y)\. 



(A6) 



Using (A6) we can estimate the off-diagonal elements 
given by (A4). Requiring these to be small compared 
to subband spacing 7i^/m£^, leads to 



VZg < 1. 



(A7) 



In summary, Eqs. (Al) and (A7) are the conditions 
for the validity of the Born-Oppenheimer approximation. 
Note that the same conditions can in fact be derived un- 
der less restrictive assumptions than Eq. (A5). 

If we assume that Zg is statistically the same as the 
surface imaged by atomic force microscopy, Eq. (A7) 
may be rewritten 



AFM 



L2 



< 1, 



iAFM 



< 1. 



(A8) 



For the control sample studied by Anderson et al. ||] 
Aafm = 1-35 A, L = 0.055 /xm. A/ = 100 A; for the 
textured sample, Aafm = 7.5 A, L = 0.09 /xm. A/ = 
70 A. Taking ^ « 30 A, it is easily seen that the Born- 
Oppenheimer approximation is applicable. 



APPENDIX B: COOPERON PATH INTEGRAL 

The purpose of this Appendix is to calculate the inte- 
gral, Eq. (36). It is useful to first calculate 

/ + 00 
dyi...%„i^(")(2/,yi)...Z?(")(y„,y'); 
-oo 

i^(o) (y, y') ^ J_e-iy-v'?-a{y+y'? . (bi) 



(A4) It is easy to show by induction that Z)'") is of the form 



1 1 



D(")(y,y') = J^^e-<^r.[y~y') -br.(y+v') (33) 
with a„&n = g. By completing squares we find 



Z?(«+l)(y,y') 



+ 00 



dyrD^°\y,y^)D^^\y^,y') 



1 1 



provided a„6„ = g. Here 

1 1 1 



(B3) 



a„+i an + g bn + 1 
1 1 1 



bn+l an + 1 bn + g ' 



An+l = ^n\/an + 6„ + 1 + 5. 



(B4) 



Using Eq. (B4) we can verify that a„+i6„+i = g ensuring 
£)("+2) g^igQ |-,g of the form given in Eq. (B2). In 
principle the recurrence relations, Eq. (B4), determine 
a„, bn and An', but in practice it is easier to follow a 
different method. 

Experience with paths integrals suggests that the 
eigenfunctions of the kernel D'"' are Hm{V^y)e~"^ 
with a suitably adjusted — of the same form as the eigen- 
functions of a harmonic oscillator. To verify this conjec- 
ture and to fix a we use a generating function for the 
Hermite polynomials 

00 jn 

exp (-a;2 + 2xV2^y - ay^^ = ^ ^iJ„(%/2^y)e-"^' . 

(B5) 



Upon completing squares we find 

J dy'D'-"^ {y, y') exp(-x2 -f 2xV2^y' - ay'^) 



i/„(2gi/4y)exp(-2V?y") (B6) 
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provided a = 2^^. This condition is imposed by re- 
quiring that the coefficient of the term in the second 
Une of Eq. (B6) should be a. Comparing Eqs. (B5) and 

(B6) we conclude that iJ,„(2,gi/-*y)e-2yiy' 

are eigenfunc- 

tions of £>(°) with eigenvalue (1 - v^)"/(l + v^)™+^ for 
m = 0,l,2,... 

Now let us evaluate 

/+00 
dy'D^^\y,y') exp(-2Vffy'' - ar^ + 4g'/^xy) 
-oo 



and = B±/Be. In the Landau gauge iJ^"^ depends only 
on (x — x') and has the Fourier transform 



1 



/Trar, 



1 + V9, 



xy 



1 ^/^^a^ {an - 



E 



m=0 



m! An {an + s/g)"'+^ 



x^H^{2g^l^y)e-^^^y' 



E 



(1 



/^\m(rt+l) 



The integral in the first line of Eq. (B7) has been ana- 
lyzed in two ways. First we use the ansatz for Eq. 
(B2), perform the integral by completing squares, and 
expand the result using the generating formula for the 
Hermite polynomials, Eq. (B5). The results are given in 
the second and third lines. Alternatively, we expand the 
exponential in the first line in terms of the cigcnfunctions 
of using the generating formula, Eq. (B5), and note 
that the eigenfunctions of are also eigenfunctions 
of D^"^ but with the eigenvalue raised to the [n + 1)"' 
power (considered as an integral operator, amounts 
to n + 1 repeated applications of D^^^). The result is 
given in the fourth line. Comparison of the third and 
fourth lines reveals 

-(l + Vg)"+i + (l-Vg)"+i- 



K = 



V9 



(1 + - (1 - 

(1 + - (1 - Vg)"+i 

(1 + + (1 - 



n+l 



X [(i + v5r+'-(i 



Nn+l] 1/2 



(B8) 



v^r+i-(i-V5r 

These expressions are the solution to the recurrence re- 
lation in Eq. (B4) with the initial condition ao = 1, 6o = 
9-,Ao = 1. 

In summary the integral D^"^ defined in Eq. (Bl) is 
given by Eqs. (B2) and (B8). 

Now we evaluate Eq. (36) for i^^(r, r'). Rescaling 
the co-ordinates we obtain 



211 



(B9) 



with 



E^^\v,y')= J dri...c/r„i;(°)(r,ri)...E(o)(r„,r'), 
^(°)(r,r') = lcxp[-\r-r'f-ip{x-x'){y + y')] (BIO) 



exp [-(a; - x'f - 0{y + y'){x - x')] 



/■+°° dk f 



^^exp.-M|±^]e^^(-') 



(Bll) 



By virtue of translational invariance 



£;(")(r,r') 



+~ dk , , 1 

— / dyi... dyn 



■n+l 



^ik{x—x') 



-[k+P{y+yi)]V4 



-lk+f}{y^+y')]V4: 



xe 



-{y-yi) 



(B7) 



-n+l 



(B12) 



^ik{x—x') 



If we shift yi ^ yi — k/2f3 

^ dyi...dy„ 

Xexp -(^{y + y^f + ^{y^ + y'f 

^(.-{y-yif+--{yn-y'r (b13) 

Here y = y — fc/2/3, y' = y' — k/20 and the fc-integral 
must be performed after the y-integrals. Evidently 



E(")(r,r') = ge^Mx-x')^(n) (^y _ k_ y, 



A 

(B14) 



£)(") is given by Eqs. (B2) and (B8) with g 13^/4. 
Performing the remaining Gaussian integral over k yields 

1 12 II 

_£;(")(!• = g-a„|r-r | -ifiiy+y )(x-x ). 

[l + (/3/2)]"+i + [l-(/?/2)]"+i_ 



an = 



[l + (/3/2)] 



n+l 



[1 - (/3/2)]" 



+1 ■ 



1 

/3 



n+l 



/3 



n+l 



(B15) 



Eqs. (B9), (BIO) and (B15) constitute the final expres- 
sion for (r, r'). 

The calculation of is similar to the solution of 

Schrodinger's equation for an electron in a magnetic field 
because the slowly varying eigenfunctions of k'^^^ obey 
the differential Eq. (33). However the derivation of Eq. 
(B15) goes beyond solution of the differential equation 
because we also obtain the short distance behavior of 



K 



K 



(n) 

Setting r = r' in Eq. (B15) we obtain 

(n-l)/. _n 1 -B+ 



\r,r') 



27r|2 B. 



1 



25, 



25« 



2(/>s, 



sinh 



nBj_ 
2^7 



(B16) 
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provided that B±/Be is sufficiently small that not only is 
B±/Be <C 1 but also {B±/Be)'^n < 1. The approximate 
expression in Eq. (B16) is the same as Eq. (41). For 
a given large Be/B± Eq. (41) is accurate only for n ^ 
{Be/B±)'^; but since K^j^^ (r, r') is exponentially small for 
n ^ Be/B±, no significant error is made by taking Eq. 
(41) to apply for all n. 



APPENDIX C: DIRECTED AREA 
DISTRIBUTION 

First let us make precise the meaning of directed area. 
Imagine an n-sided polygon drawn on a plane. The poly- 
gon is endowed with an orientation by assigning each 
edge a direction such that each vertex has one incoming 
edge and one outgoing. Roughly, the two orientations of 
a given polygon correspond to the two senses in which 
its perimeter may be circumnavigated. 

Now suppose a magnetic field is applied perpendicular 
to the plane. Roughly, the directed area is the quantity 
by which the magnetic field must be multiplied in or- 
der to obtain the flux through the oriented polygon (in 
magnitude and sign). For a simple polygon with no self- 
intersection the directed area is equal in magnitude to 
the area; the sign is positive or negative depending on 
whether the orientation of the polygon is anti-clockwise 
or clockwise viewed from the direction in which the mag- 
netic field points. For a polygon with self-intersection, 
the directed area is obtained by weighting the area of 
each infinitesimal area element by its winding number 
before adding them. Here the winding number of an area 
element is computed by drawing a vector from the ele- 
ment to the boundary of the polygon and counting the 
number of anticlockwise turns made by the vector as its 
tip circumnavigates the perimeter (clockwise turns count 
as negative turns). For a polygon with successive vertices 
at r, ri, . . . , r„_i, r immersed in a magnetic field pointing 
in the ^-direction, the directed area is given by 

a = - ^)iyi + y)- - xi){y2 + yi) + ... 



P{a) 



da /na\ 



^{x- Xn){y + yn). 



(CI) 



Let the probability density of drawing an n-sided ori- 
ented polygon with successive vertices at r, ri , . . . , r„_i , r 
be (n/7r"-i)exp[-|r-ri|2-f ...-f|r„_i-r|2]. Eqs. (CI) 
and (BIO) show that 



1 

^("-i)(r,r) = — / daP(a)exp(i2/3a). (C2) 

Here P{a) = probability distribution of the directed area 
for the Gaussian polygon ensemble defined above. By 
inverting the Fourier transform in Eq. (C2) and making 
use of Eq. (B15) for (r, r') we obtain 



oo 

+00 



27r V 2 / 



TT 2 I 27ra 
— — sech 



(C3) 



The approximate form of the integrand in the second line 
above is valid only for a <C 1/ ^/n•, but it may be extended 
to all a since the integrand is already exponentially small 
for a 3> 1/n. 

Eq. (C3) is the formula for the directed area distribu- 
tion. It coincides with Eq. (41) if we take the step size 
of the random walker to be \/2le instead of 1 as we have 
done in this Appendix. 

The distribution of directed area for Brownian motion 
was analysed by Levy | |3C| ] and restudied more explicitly 
in connection with electron transport by Argaman et al. 
pit . Here however we require the area distribution P„ (a) 
for a random walk with a finite number of steps n. The 
large n behaviour of (a) was recently obtained by Bel- 
lissard et al. for random polygons on a square lattice by 
an interesting application of non-commutative geometry 
P2[. Their result agrees with the large n limit of eq (C3). 
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